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Abstract 
>v 

We develop a general theory for active viscoelastic materials made of polar filaments. This theory 



is motivated by the dynamics of the cytoskeleton. The continuous consumption of a fuel generates 
a non equilibrium state characterized by the generation of flows and stresses. Our theory can be 
applied to experiments in which cytoskeletal patterns are set in motion by active processes such as 
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I. INTRODUCTION 

Molecular biology has scored an impressive number of success stories over the last forty 
years. It has provided us with a vast knowledge about molecules at work in living systems, 
about their involvement in specific biological functionalities and about the structure of their 
chemical networks [l|. Yet one is still unable to take advantage of this knowledge for con- 
structing a comprehensive description of cell behavior. Assuming that we had all desired 
molecular informations, the computer time required to describe meaningful cell behavior 
would be totally prohibitive. Meanwhile there is a clear need for a generic description of 
cells. 

An alternative approach consists in identifying a reduced number of key characteristics 
and construct a phenomenological description of a "simplified" but relevant cell. Again 
this task is currently too complex at the scale of a whole cell, but it can be envisioned for 
some of its constituents. A good example is provided by the description of membranes: 
for length scales of the order of a few tens of nanometers up, it can be considered as a 
fluctuating surface on which a few densities are distributed and through which a few other 
densities p_ passively 0r actively QQ. The faction of the theory has been g oing 
on for more than thirty years. It sheds light on the physics of membrane shape, topology 
changes such as budding, equilibrium and non equilibrium fluctuations, membrane adhesion, 
long range interactions on membranes etc... It is still under development to include out of 
equilibrium features such as lipid and protein exchange with the bulk, but much has been 
learned already. 

Efforts towards understanding the cytoskeleton are more recent. They have been focused 
on the description of its passive visco-elastic properties which are now fairly well understood 
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Such materials which can 



in terms of gels made of cross-linked semi-flexible polymers 
be prepared in vitro, are equilibrium systems which obey conventional thermodynamics. In 
eukaryotic cells, the problem is however qualitatively new: the cross-links can be made of 
motor proteins which have their own dynamics driven by chemical energy. Experiments, 
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rave shown that such systems have a rich and corn- 
One can grasp the degree 



of complexity with the remark that cross-links define distances, in other words they define 
a metric. The cross-linking agents being motor proteins they move and the metric evolves 



with time. 

In fact the problem is even more complex since cytoskeletal filaments are polar and out 
of equilibrium: they polymerize at one end while depolymerizing at the other end. Such 
a process called treadmilling is well known to biologists. The description of eukaryotic cy- 
toskeletal gels should include all these features. We call these gels and more generally all 
gels working in the presence of a permanent energy consumption "active" gels. On large 
length scales and long time scales, the properties of complex materials can be captured by a 
generalized hydrodynamic theory based on conservation laws and symmetry considerations. 
Such theories have been very successful in describing complex fluids such as liquid crystals, 
polymers, conventional gels, superfluids etc. J2l|, UA, l23j |. For instance all liquid crystal 
display devices can be described by such theories. More recently active liquid crystals have 
been considered within this same logic J24J. Here, we develop a hydrodynamic theory of 
active gels. An example of such gel is the network of actin cytoskeletal filaments in the 
presence of chemical processes which locally induce filament sliding and generate motion. 
These active processes are mediated by motor proteins which hydrolyze a fuel, Adenosinet- 
riphosphate (ATP). Since cytoskeletal filaments are structurally polar, each filament defines 
a vector. This filamental structure can on large scales give rise to a polarity of the material 
if filaments are aligned on average. 

We develop the hydrodynamic theory of active polar gels systematically in several steps, 
following standard procedures. First, we identify the relevant fields and write conservation 
laws for conserved quantities. We identify the generalized fluxes and conjugate forces in 
the system. These fluxes and forces define the rates of entropy production and dissipation. 
Using the signature of forces and fluxes with respect to time reversal, we define dissipative 
and reactive fluxes. These fluxes can be expressed by a general expansion in terms of forces, 
by writing all terms allowed by symmetry and by respecting the signatures with respect to 
time reversal. By keeping the relevant terms to lowest order, this finally results in generic 
dynamic equations which are valid in the vicinity of thermal equilibrium. We then provide 
examples by analyzing the active gel mode structure and by discussing the spontaneous 
dynamical behavior of topological singularities such as disclinations (asters and vortices) in 
two dimensions. Eventually we discuss the merits of our approach, its current limitations 
and ways to extend its domain of validity to the more relevant far from equilibrium regime. 
Some of the results discussed here have been presented in Ref. 



II. FIELDS, DENSITIES AND CONSERVATION LAWS 

We develop the hydro dynamic description of active gels starting with conservation laws. 
The number density of monomers in the gel is denoted by p. The gel is created by poly- 
merization of filaments from monomers. The monomers have a density p^ in solution. The 
treadmilling process and the active stresses in the gel create a flow of the actin gel monomers 
with a local velocity v. The gel current is then convective and mass conservation implies 

that 

dp 

-£ + V(vp) = -k dP +k p 5 s (1) 

Here k d is the depolymerization rate; we assume, as seems to be the case for actin, that the 
depolymerization which occurs mostly at the branching points of the gel, has a rate propor- 
tional to the local density. In some cases treated below, we will consider for simplicity that 
depolymerization does not occur in the bulk but at the surface of the gel. In many situations, 
filament polymerization is highly controlled and localized for example at the surface of the 
gel. This is taken into account in Eq. ((TJ by introducing the surface polymerization rate 
k p and 5s denotes a Dirac-like distribution which is non- vanishing only at the gel surface S. 
Similarly, we write a conservation law for free diffusing monomers, 

J^_ + V ^=k d p-k p 5 s (2) 

here, we have introduced the diffusive current j^ of free monomers. 

Active processes are mediated by molecular motors. The concentration of motors bound 
to the gel c^ is an important quantity to characterize the effects of active processes in 
the gel. Assuming for simplicity that the total number of motors is conserved, we write 
conservation laws for c^ and the concentration of freely diffusing motors in the solvent c^ 
which read 

Q c (m) 



dt 

d c (b) 



+ vj (m) = k oS c^-k onP {c^y 



<)t +Vc«v + Vj« = -k oS c^ + k on p(c^) n (3) 

The attachment and detachment rates of motors to and from the gel are characterized by 
the chemical rates k on and k g. Here, we have taken into account that bound motors are 
convected with the gel. The current of free motors is j^ and we denote j*- 6 -* the current of 
bound motors relative to gel motion. In general, binding of motors to the gel is cooperative, 



and cannot be described as a second order reaction: groups of motors could bind together. 
We use here an n th order chemical kinetics where the rates is proportional to (c^) n . 

A final important conservation law is momentum conservation. In biological gels on scales 
of micrometers, inertial forces are negligible and momentum conservation is replaced by the 
force balance condition which reads 

9 Q (^-nM + /r = o (4) 

where /| xt is an external force density. Locally, there are two forces acting on the gel, the 
total stress tensor a^a and the pressure II. 

The dynamics of the system is specified if the flow velocity v and the currents j^ a \ j*" 1 -* and 
j( fe ) are known. The physical description of the currents is discussed in the following sections. 
Furthermore, we have to take into account the polar nature of the gel. Individual filaments 
are rod-like objects with two different ends which therefore have a vectorial symmetry. If 
the filaments in the gel are on average aligned, the material is oriented. We introduce a 
polarization field p to describe this orientation. The field p is defined by associating with 
each filament a unit vector pointing to one end. The vector p is given by the local average 
of a large number of these unit vectors. 

III. CONSTITUTIVE EQUATIONS 

The fluxes of monomers and motor molecules are generated by forces which act on the 
active gel and induce motion. In this section, we identify the relevant forces and derive 
general flux-force relations. These relations define the material properties, and characterize 
how the system reacts to different types of generalized forces. Of particular significance for 
our theory is the existence of active processes mediated by molecular motors. In general, a 
chemical fuel, such as Adenosine-triphosphate (ATP), is used as an energy source. Motor 
molecules consume ATP by catalyzing the hydrolysis to Adenosinediphosphate (ADP) and 
an inorganic phosphate and transduce the free energy of this reaction to generate forces 
and motion along the filaments. The energy of ATP is also used for polymerization and 
depolymerization of the filaments. The presence of the fuel is equivalent to a chemical 
"force" acting on the system. We characterize this generalized force by the chemical potential 
difference A/x of ATP and its hydrolysis products, ADP and phosphate. When A/x vanishes 



the chemical reaction is at equilibrium and there is no energy production. When A/j, > 0, a 
free energy A/x is consumed per hydrolyzed ATP molecule. 

Constitutive equations are obtained by first identifying the fluxes and the corresponding 
conjugate generalized forces and then writing a general expansion for the fluxes in terms of 
the forces. We study in this paper an active gel close to thermodynamic equilibrium and 
limit the expansion of fluxes in terms of forces to linear order as in a standard Onsager 
theory. We thus describe the linear response of the gel to generalized forces. We do this in 
the most general way and write in the flux-force expansion all terms which are consistent 
with the symmetries of the system. 

A. Fluxes and forces 

We first discuss the rate of entropy production in the active gel. A change in the free 
energy F = U — TS per unit time can be written as 

F = -Jdh {a^d a v p + h a p a + A/xr + dj® /x« } (5) 

where the "dot" denotes a time derivative. The total deviatory stress tensor a^ is in general 
is not symmetric; it is conjugate to the velocity gradient d a vp. The field conjugate to the 
order parameter p a is the functional derivative of the free energy F of the gel at thermal 
equilibrium, h a = —j^~, where the functional derivative is taken for constant deformation, 
temperature and number of particles. The current conjugate to the field h a is the convected 
time derivative of the polarization p a = -Jjp + vpdpp a . The chemical force A/x is conjugate 
to the ATP consumption rate r which determines the number of ATP molecules hydrolyzed 
per unit time and per unit volume. Finally, jj> a \ ^ m ' and /J® are the chemical potentials 
of free monomers, free motors and of motors bound to the gel, respectively. 

Eq. (J5J) does not take into account the translational and rotational invariance of the active 
gel, and the variables used are therefore not the proper conjugate fluxes and forces. Indeed, 
as shown in reference J32[ , since the free energy does not change under pure translations and 
rotations of the gel, ignoring surface terms, we can rewrite Eq. (|5J) as 

F = -Jdh {a a0 u a0 + h a P a + Afir - j«d a /x W } (6) 

The total stress a tot has been decomposed into a symmetric part a a p (with a a p = crp a ) and 
and anti-symmetric part which is due to the torque exerted by the field h a on the order 



parameter p a : 

cS? = ^a/3 + 2 (P* h ~ Pf* h <*) ( 7 ) 

Here, u a p = l(d a vp + dpv a ) is the symmetric part of the velocity gradient tensor. The 
current conjugate to h a includes a rotational contribution coming from the antisymmetric 
part of the stress tensor, it reads 

D dv 

Pa = JZ-Pa = -KT + (Vjd^Pa + U a pPp (8) 

where we have used the corotational time derivative of the vector p a , ui a p = \{d a vp — dpv a ) 
being the vorticity tensor of the flow. 

We can read off Eq. © the following pairs of conjugate fluxes and forces: 

flux <-» force 



(9) 
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The rate of change of the free energy can be divided into reversible and irreversible parts 
F = F Tev + F irr where F rev = U and F irr = — TS. We therefore decompose the fluxes into 
a reactive part and a dissipative part. They are characterized by their different signatures 
with respect to time- reversal. Note, that the generalized forces have well-defined signatures 
with respect to time reversal: the rate of strain u a p is odd under time reversal, while h a and 
A/x are even. We write 
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(10) 

The dissipative fluxes have the same signature under time reversal as their conjugate forces, 
while reactive fluxes have the opposite signature under time reversal. Therefore, a%,8 * s °dd 
under time reversal, while r d and P d are even. The reactive parts have correspondingly the 
opposite signatures. As we show below, the currents j® do not have reactive parts. With 



this decomposition, the rate of entropy production TS = — F irr is given by 

TS = J dh {r d Afi + P d Ji a + otpUafi - J^d^ } . (11) 



[ T dt(r r A f i + P r a h a + a r af3 u a p) = . (12) 

•J U 



For a changing state of the system which is periodic in time with period T, the internal 
energy is the same after one period. Therefore, in this case 
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These relations follow from the signatures of dissipative and reactive currents under time 
reversal. 

B. Maxwell model 

In the absence of a net gel polarization p and for a passive system with A/z = 0, we 
describe the viscoelastic gel by a Maxwell model. In order to keep our expressions simple and 
to focus on the essential physics, we assume that the ratio of the bulk and shear viscosities 
is 2/d where d is the space dimension. The Maxwell model relates the stress tensor to the 
strain rate according to 

f 1 + T— J <J aj 3 = 27]U af3 (13) 

Here, rj denotes the shear viscosity and r the viscoelastic relaxation time, and E = rj/r is 
the elastic modulus of the gel at short times. In Eq. I)13j). the time derivative of the stress 
tensor taken in a reference frame moving with the material flow must be used. For a tensor, 
this implies that in the laboratory frame convective terms due to the fluid flow as well as 
the rotation of the reference system due to the vorticity of the flow need to be taken into 
account. D/Dt denotes the corotational derivative of a tensor given by 



Da al3 _ da a/ 3 



+ (w 7 (9 7 )cr Q/3 + [w a7 cr 7/3 + ^/3 7 cr 7a ] (14) 



Dt ' dt 

Note that we include the geometrical non-linearities and that we use the most general ver- 
sions of what is called a "Convected Maxwell Model" ||23|. For the Maxwell model given by 
Eq. (fI5|L we can specify the reactive and the dissipative contributions to the stress tensor. 

(i-^-RsKu = ^nuafi (is) 



Dt 2 ' 



D 



Otf = - t dI°U ( 16 ) 



With these relations, we have a a/ 3 = a^p + a a/3 an d the dissipative and reactive parts of the 
stress differ in their signature with respect to time reversal. 



C. Dissipative fluxes 

We now generalize the relations (|15|) and (J16JI expressing dissipative and reactive fluxes 
in terms of generalized forces for an active polar gel. First, we write expressions for the 
dissipative fluxes. To linear order we find 

( 1 ~ r2 ^) Cr ^ = 2r]Ual3 

- T J^^(P ah P + p P ha "> + "iPi h i ( 17 ) 

D . _ 

+ T—{h>iP l3 U a p + h> l p a Ufl l 3) (18) 

r d = AAfi + \ lPa h a + \p a d a fi (b) . (19) 



Note that the first expression is a generalization for Eq. ()15j) of the Maxwell model. For a 
sake of simplicity, we ignore the anisotropy of the viscosity and consider that all translational 
viscosities are equal. This anisotropy could be introduced as in the classical description of 
the hydrodynamics of liquid crystals. Couplings to A/x cannot appear in the equation of 
the dissipative stress because it transforms differently under time reversal than u a p. A 
term coupling a d a to the time derivative of h a can appear. It is required by Onsager 
symmetry relations since a corresponding coupling term with coefficient V\ does occur in the 
equation for P d as shown in Appendix B. The dissipative coefficients 71 and A characterize 
the coupling of P d and r d to h a and A/i, respectively. Here, 71 is the rotational viscosity 
which appears in classical nemato-hydrodynamics. Because of its signature with respect to 
time reversal, cross-coupling terms involving u a p do not appear in the expression r d . The 
expression for P d contains a terms coupling P d to the time derivative of u a/ 3 as derived in 
Appendix B. Cross-coupling terms occur also in the expressions for P d and r d . Since P a is a 
vector and r a scalar, we need a vector in the system to couple these equations. If the system 
is polar, this vector is provided by p. Therefore, the cross-coupling term characterized by Ai 
involves p a . Note that because of the symmetry of the Onsager coefficients, the coefficients 
Ai V\ and V\ are the same in different equations. The last term in the expression of r d is a 



crossed term with the current of bound motors discussed below. This term is also imposed 
by the Onsager symmetry relations. We still have to specify the currents j™ , j^ 1 ' and jj® 
of monomers and motor molecules which do not have reactive components: 

j>) = -D^d a p {a) (20) 

j(m) = _ D (m) daC (m) ^ 

£) = -D®d a c® + \ Pa Afi (22) 



The first two expressions are simple diffusive currents with diffusion coefficients D^\ where 
we have assumed for simplicity that couplings to h a can be neglected and that free diffusion 
is not influenced by A/i. The expression (I22|) describes bound motors which are convected 
with the gel. In addition to being convected with velocity v, motors move in a direction 
given by the filament orientation if A/i is nonzero. This directed motion is characterized 
by the coefficient A; the motor velocity v m on the filaments is such that c^'v™ = Xp a Afi. 
Finally, velocity fluctuations of motor motion are captured by a diffusive term with diffusion 
constant D^ b \ 

D. Reactive fluxes 

Reactive fluxes are expressed in terms of the forces following their signature under time 
translation invariance. Generalizing Eq.([16|) for the reactive stresses of the Maxwell model, 
we note that here cross coupling terms with the forces h a and A/x are possible. Taking into 
account the tensorial structure of the stress tensor, we write 






D <P +A 

+ A af3 



- (A[ip a pp - (An8, 



a/3 



Dt 

+ -^(Pahp + Pph a ) + VlP~fh^5 a p , (23) 

where we have introduced the phenomenological coefficients (, (, (', v x and v\. The tensor 

A al 3 = v 2 (u ai a 1 p + a ai u ll3 ) + vzU 11 a al3 + v A u 11 a 56 8 a p 

+ u 5 a^u a + f&u^af^ap . (24) 

contains nonlinear reactive terms to lowest order, resulting from the geometry of the flow 
field with corresponding phenomenological coefficients i^. Similar coefficient have been in- 
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troduced in the so-called "eight constant Oldroyd model" in rheology [26j. The term pro- 
portional to r on the right hand side of Eq. (J2"3j) assures compatibility with the Maxwell 
model in the absence of polarization and chemical fuel. For the reactive parts of P a and r 

we write 

D 2 

(1 - r2 £^2) P a = -WapPP ~ U l U WVa (25) 

r T = (PaPl3U a p + (U aa + ('paPaUp/3 (26) 

Here we have written cross coupling terms of P r with the rate of strain u a @. The linear 
response matrix of reactive terms is antisymmetric, therefore the same coefficients as in Eq. 
(J23j) appear, however with opposite sign. The remaining Eq. (J26j) is constructed in the same 
way with cross coupling coefficients that have been introduced in Eq. (|23|) . No reactive 
cross-terms between Eqns. P3j) and (J2B1) exist. 

E. Dynamic equations 

Using the expressions for dissipative and reactive fluxed discussed above, we now write 
general hydrodynamic equations for the active viscoelastic and polar gel. Adding the dissi- 
pative and reactive parts, we find 

2r]U a/3 = (l + T—j {a a p + (A/ipaPp 

- -^■(jPcchp + pi 3 h, a )-p 1 p 1 h, 1 8 0l p (27) 

-(l-T—)(viU a 0P/3 + l>iUi3i3Pa) (28) 

r = (PaPf3Uaf3 + (U aa + C'PaPcxUpp 

+ AAfi + \ lPa h a . (29) 

In these equations, we have included geometric nonlinearities but have restricted other terms 
to linear order for simplicity. Also, we have neglected chiral terms which in principle exist in 
cytoskeletal systems but which are expected to be small. These equations are complemented 
by the force balance condition (0) . 
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Eq. (|27j) generalizes the expression of the stress tensor in the Maxwell model to active 
systems with polarity. Indeed, even in the absence of stresses, the active terms proportional 
to A/x generate a finite strain rate. Similarly, if all flows are suppressed, the active terms 
generate a nonzero stress tensor. Thus, the hydrolysis of ATP can generate forces and 
material flow in the gel via the action of active elements such as motors. These terms are 
characterized by the coefficients (, (' and ( Similarly, we find active terms in the polar- 
ization dynamics given by Eq. (f^Hj) described by the coefficient Ai. Furthermore, material 
flow couples to the polarization dynamics via the coefficients v\ and V\. The rate of ATP 
consumption r is primarily driven by A/i and characterized by A. However, it is also coupled 
to the fluid flow and to the field h acting on p. Note, that in addition surface terms can 
be important. For example, if filaments polymerize at the gel surface with a rate k p (see 
Eq. (JU), there is an additional contribution rs to ATP consumption which is localized at 
the surface. In the following, we discuss situation where surface effects can be captured by 
introducing appropriate boundary conditions. 

IV. HYDRODYNAMIC MODES 

We now calculate the linear relaxation modes of an active polar gel. They can be obtained 
from the linear response of the change in the gel density to externally applied forces. We 
focus on the situation of a fully polarized gel far from an isotropic polar transition with the 
polarization vector lying in the x-direction, i.e., p = e x is constant. We assume that the 
gel is treadmilling, which implies that it moves at constant velocity v = Voe x , while the gel 
mass is conserved. The treadmilling is due to polymerization and depolymerization of the 
filaments at the gel surfaces located at x = x + and x = x~, respectively. In the following, 
we consider the case of an infinite system where x + are located at +00 and x~ to — 00 while 
the treadmilling velocity vq remains constant. 

We restrict ourselves to modes for which the fields vary along the x-direction which is 
the treadmilling direction, while the system is assumed to be homogeneous in the other 
directions. In this case the dynamic equations ©-© for the densities of polymerized actin, 
bound and free myosin become effectively one-dimensional and read 

d tP + d x (pv) = (30) 

d t c {m) - D8 2 x c {m) = k oS c {b) - k onP {c {m) ) n (31) 
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d t c^+d x (c^v + \Ap) = -hose® + k oaP (cW) n . (32) 

The actin monomers which are not part of the gel diffuse across the system and control by 
their concentrations the treadmilling rate. The mechanical constitutive equation (|27j) takes 
the form 

2r]d x v = [1 + rd t + rvd x ] {a + ( A/x) (33) 

whereas the force balance in the presence of an external stress a ext is expressed by 

d x (a + <x cxt ) = d x Tl . (34) 

In these equations, the coefficients A and (, which are the Onsager coefficients introduced 

in the previous section, characterize the directed motion of motors along the filaments and 

the actively generated stress in the gel (we have replaced here ( + £' + ( by £.). In a 

situation far from equilibrium, nonlinearities become important. We therefore linearize 

around the homogeneous steady state of the above equations given by p = p , c^ = c , 

C M — Cq , v = v , and a = cr , where k o ffC = /c on Po( c o ) n — Po v o an d °"o — — C^/x. We 

consider a small perturbation of this state, i.e., p = po + Sp and analogously for the remaining 

quantities. To leading order, the active coefficients can be expanded as ( = (o + ( p 5p + ( c 8c^ 

and A = Ao + \ p 5p + \ c Sc^ to implement such nonlinearities which become relevant in all 

realistic situations. Here, we have neglected the diffusion of bound motors and the diffusion 

constant D^ has been set to zero. 

The linearized dynamical equations can be solved by Fourier-transformation in space and 

time. The mode characterizing the variations of the density of unbound motors, c^ (q, to) as 

a function of spatial wave vector q and temporal frequency u, can be obtained from Eq. (p?T|) 

and is given by 

(m) ( N _ fc off c(*°(g,a;) ~ fcon(4 m) ) n p(g,^) , . 

too + Dq z + l c 

where T^ 1 = nk on po(cQ ') n ~ 1 is a chemical relaxation time. Using this result as well as the 

linearized conservation equation 

iup(q, uj) + iqp v(q, uj) + iqv p(q, uj) = (36) 

we obtain for the distribution of bound motors 

c{b) M^-^pJmM (37) 

A(q,uj) + iqX c po 
13 



where \ p = \ p App Q , \ c = A c A/i, and 



A(q, u)=iu} + iqv + - fc off 

The mechanical equations ffTTT^f) and (J33|) then lead to 



(38) 



a ext (g,w) = 2/i(g,^ 



Po 



(39) 



where the modulus p(q,u>) characterizes the effective material properties in response to 
density variations. It can be determined from the Eq. 



(1 + it(uj + qv )) [2//(g, u>) - XPo] = 2iEr(u + qv ) 

A(q,u)c -iq\ p 



+ (1 + iuoT + iqv r) 



C + Cc 



A(q,u) + iq\ c 



(40) 



Here, we have used the equation of state and have defined the inverse compressibility x — |r- 
Furthermore, we have introduced ( p = ( p Afip as well as ( c = ( c A[i. 

The relaxation modes of the system can be identified by the (complex) values of the 
frequency to n , for which the modulus vanishes /j,(q,u n ) = 0. Indeed, the density relaxes as 
p(q) ~ e luJnt and instabilities in the system occur if the real part of iu n becomes positive. 
For our present system, we find that three independent relaxation modes exist, with n = 1 
to 3. Up to leading order in q, the corresponding dispersion relations are given by 



lUOi 

iuj 3 



icq + dq 2 

-T~ l - k oS 

1 XPo + CcCo + C P 



t2E + xpo + CcCq + ( P 
The first mode is propagating with velocity c given by 



(41) 
(42) 
(43) 



1 r K, 



off 



v + 



(XPo + Cp)A c - C c A f 



(44) 



XPo + (p + CcCo 

Note, that the propagation velocity c is proportional to the fraction of bound motors (1 + 
T c k g)~ x . The value of the coefficient d can be either positive or negative, depending on 
the value of the parameters. If it is positive, the gel is unstable and self-organizes into 



propagating density profi 



of active polar fibers 
in several cell types 



es. Such solitary waves have been found in theoretical descriptions 
|l6(. Experimentally, actin waves have been recently reported 



JlaJla- ] 

y, baboj 



In the present context, there is only one propagating 
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mode, because space inversion symmetry has been broken due to treadmilling and to the 
polarization of the gel. The polarization couples to the filament and motor densities only 
through the active current of bound motors. If instead of motors, we have passive cross- 
linkers, i.e., A = and in the absence of treadmilling, i.e., Vq = 0, space inversion symmetry 
is restored. In this case c = and 

d = ~D-^— , (45) 



T" 1 + k, 



- c 



off 



so that the mode is no longer propagating, but diffusive and always stable. Note, that 
d is proportional to the fraction of unbound motors T c k s/(1 + T c k ff). The existence of 
propagating waves even in a viscous environment seems to be a general feature of active 
media and has already been reported in an other context by Ramaswamy and coworkers 
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3l|. 



The second mode is a chemical relaxation mode describing the binding and unbinding 
of molecular motors to the polymerized gel, whereas the third mode describes the stress 
relaxation towards its stationary value. 

V. DYNAMIC POINT DEFECT IN TWO DIMENSIONS 

As an example of active behavior in two dimensions described by Eqns. ()29|1 . we consider 
in this section the dynamics of point defects in the vector field p. First, we consider the 
passive equilibrium state with A/x = where all fluxes vanish. Subsequently, we determine 
stationary active solutions for finite A/i and determine their stability. We obtain a complete 
diagram of states for asters, vortices and rotating spirals. 

A. Asters and vortices at thermodynamic equilibrium 

In order to study point defects in two dimension, we consider for simplicity the situation 
where the orientation of p varies but the modulus is constant and we impose p 2 = 1. In 
this case, the free energy is given by the standard expression for a polar liquid crystal [32]: 



d x 



y(V- P ) 2 + (p-Vp) 2 



+ fcV-p /lyp 2 



(46) 
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where K\ — K and K% — K + 5K are the splay and bend elastic moduli and we have 
introduced a Lagrange multiplier h\\ to impose the constraint p 2 = 1. The coefficient k 
describes the spontaneous splay allowed by symmetry for vector order. Note that there is 
no twist term in two dimensions. 

For A/i = 0, the system reaches an equilibrium steady state where all fluxes and forces 
vanish: a a p = 0, u a p = 0, P a = 0, h a = and r = 0. The equilibrium orientation of the 
polarization p, corresponds to a vanishing orientational field h a = —5F/5p a = 0. In order 
to describe point defects, it is convenient to introduce polar coordinates (r, 8) and the angle 
if)(r,8) which characterizes the components of the vector p: (p r = cosif), pe = sinip), see 
Fig. 1(a). The components of the orientational field h a in cylindrical coordinates can be 
written as h r = h\\ cosip — h± sinif), and hg = h\\ sin^ + h± cosi/j, where we have introduced 
the components h\\ and h± parallel and perpendicular to the direction of p. 

Considering for simplicity rotationally symmetric fields describing point defects with if) = 
ijj(r), we obtain 

2 



r , \k (i d V 

Z7r / dr r { — — — r cos w 
J )2\rdr r J 



K + 5K (I d . , . 

where we have ignored the spontaneous splay which leads to a boundary term. 

The perpendicular field is obtained by taking the functional derivative of the free energy 
h ± = -5F/5ip 



h ± = (K + 5K cos 2 ip) 



r 



5K 

sin 2w 

2 Y 



\i + ^' 2 



(48) 



where the prime indicates derivation with respect to r. The value of the Lagrange multiplier 
h\\ (the longitudinal field) is chosen in such a way that the condition p 2 = 1 is satisfied. 

Four types of topological defects of charge one are possible, see Fig. 1(a) and (b). If 
we assume that boundary conditions are always chosen to allow for solutions with constant 
if), these solutions corresponds to asters with if)(r) = and if)(r) = tt and to vortices with 
if)(r) = ±7r/2. The linear stability of asters and vortices against angular perturbations Sif>(r) 
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is described by 

7i dSip 



l + 7l (z/ 1 ±l) 2 /(4r/) dt 

{K + SK cos 2 ip ) I — - - - — UV 



d 2 1 d 



x cir 2 r dr 
- <5fs: cos(2^)-^(JV(r) (49) 

where the minus sign corresponds to an aster and teh plus sign to a vortex. If all the 
eigenvalues of the linear operator on the right hand side are negative, the defect is stable, 
and otherwise it is unstable. (We may note that this operator has the form of the negative 
a quantum Hamiltonian of one particle.) Asters are stable for positive 6K while vortices are 
stable for 5K negative. For the special case 5K = 0, Eq. (}4T?|) becomes jidijj/dt = KAi/j, 
where A denotes the Laplace operator, and defects for any value ip(r) = ip are stable. These 
defects are called spirals since the direction of polarization follows spirals with an equation 
in polar coordinates given by 

r(9) = r exp[cot(ij )9} , (50) 

see Fig 1. (c). 

B. Non-equilibrium steady states 

We now discuss the effect of a point defect in two dimensions if motors in the system are 
active and A/x 7^ 0. In this situation, spiral defects start to rotate and there are nontrivial 
flow and stress profiles. We assume that the system is incompressible, i.e., u„ + uee = 0. 
Using the expression of the strain rate tensor given in appendix, this imposes that v r = a/r. 
The absence of singularity of the radial velocity for small r implies that a = 0, therefore 
v r = and u rr = u ee = 0. 

The Onsager Eq. (|2K|) for the polarization rate can then be written in cylindrical coordi- 
nates. In a steady state, P a = 0, we obtain at linear order: 

h\\ 
z^i ii r g sin 2-0 = — + AiA/i (51) 

7i 

(uicos2tp — l)u r e = — . (52) 

7i 
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where we have expressed the orientational field h a in terms of its parallel and perpendicular 
coordinates. Similarly, Eqns. ()27| ) -([29| ) can be rewritten in polar coordinates. Using u rr = 
uqq = and only taking terms to linear order in the generalized forces into account we find 
that the steady state obeys 

= a„ + CA/icos 2 ^+ (C + C')A/i 



-V\ /i||Cos 2 ?/> /ij_ sin 2?/> 

—v\h\\ 

aee + (Afi sin 2 ^ + (C + g)A(i 



(53) 



hu sin 2 ip H — h± sin 2ip 



2r]u 



r9 



—V\h 



a r e H — A// sin 2ip 



2 1 



h» sin 2ip + /ij_ cos 2^ 



(54) 



(55) 



In addition to the dynamic equations Eqns. (J2*7j) - (|2^j) . the force balance Eq. (HJ) has to be 
satisfied. The 9 component of the force balance equation ( Eq. (jA6|) of appendix EI) is solved 
by a^g = B/r 2 where B is an integration constant. Since a^ff must not diverge, B = and 
therefore a^f = and a r g = h±/2. The components a rr and aee of the stress tensor follow 
from Eqns. (J53)) and (|5"2j) . These stresses together with the radial component of the force 
balance equation (Eq. ()A5j) of appendix \K§ determine the pressure profile II required to 
ensure the incompressibility of the material. 



C. Asters, vortices and spirals 



We now determine the stationary solutions for point defects in non-equilibrium states. 
From Eq. (joTj) we find an expression for the parallel component of the orientational field 



h\\ = r yi(i'iu r esm2ip — AiA/i) 

Inserting this expression in Eq. ()55|) . we obtain the non diagonal strain rate 

( A/i sin 2ip + h± ( 1 — v\ cos 2ip) 



(56) 



u r8 



Ar] + 7 X z/ 2 sin 2ifi 



(57) 
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where C = C + ^lTiAi- This equation together with Eq. (JS2J, where h± is given by Eq. (jJHj) 
determines ^ and w r 0. Eliminating M re , we find an equation for the steady state orientation 

h±(4r) + 71(^1 +l-2ui cos 2^)) = 7l sin2^(^i cos 2^ - 1)<A// (58) 

We now discuss special solutions to this equations as well as the stability of asters and 
vortices obtained at thermal equilibrium for A/x = 0. 

1. Spiral solutions for 5K = 

We first consider the simple case where there is no anisotropy of the elastic constants: 

SK = and h± = KAip. Eq. (|52jl requires that spirals with constant ipo and h± — satisfy 

cos 2^o = 1/Vi- This selection of angle, which depends on the parameter V\ is a dynamic 

phenomenon. The same angle is selected for the orientation of nematic liquid crystals in 

shear flows. A steady state exists only if \vt\ > 1. The polarization angle takes one of the 

two values 

ip = ±- arccos — (59) 

2 V\ 

in the interval [—it, 7r]. With these values of ipo, we find 

sin 2^o - 

^ = 4, + 7l ^sin 2 2^ CA/i • (60) 

The velocity is ortho-radial, along the 9 direction; it is obtained from (jA3|) 

v e = 2r(£^dr' + v^ (61) 

where v@ is an integration constant. For a finite system with radius R and with the boundary 
condition that no motion occurs at the boundary, 

T 

v e = 2u r9 rlog— . (62) 

ti 

2. Stability of asters and vortices 

Both asters with tp Q = 0, tc and vortices ipo = ±7r/2 are solutions to Eq. (|58Jl . In order to 
understand under which conditions these solutions are stable, we perform a linear stability 
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analysis, writing if)(r) = ipo + 5ip. The time-dependence of ip in the liquid limit r — > is 

given by 

dip 
7i-oT = h± -7i(^cos2^- l)u r6 (63) 

Linearizing this equation and using Eq. ()57|) we find for asters with ip = and SK > 
that £■?/.> (r) satisfies to linear order 

- . 7l — rT7 - A . —Sip = C5tP (64) 

i + Tl (i/!- 1)7(47?) at ^ v ; 

where the linear operator is given by 

^(™# + ^"^4 («) 

Here, we have defined n 2 = SK/(K + 5iT) and A; 2 = — [4 ^ .A^i^ 1 ,^^-) - An aster solution 
becomes unstable, if the largest eigenvalue of £ vanishes. The condition C5ip = is solved 
by Bessel functions of order n, Sip ~ J n (kr). For a finite system with radius R we find with 
boundary conditions 5ip(R) = that kR = z n , where z n denotes the first positive root of 
J n (z). Therefore, the critical value for the instability of asters is given by 

A /2 n \ 2 4r7 + 7 1 (z/ 1 -l) 2 



^ = -{f) '^(A-D {K + m) (66) 

In the case of vortices with %po = ±7r/2 and 5K < 0, we find 

Jl dty _ Wr, 

l + 7l (z/ 1 + l) 2 /(4r / ) dt W 

, _ ( d 2 Id m 2 2 \ 
\dr 2 r dr r 2 J 

where m 2 = —5K/K and q 2 = — , 4 (v +iv\k ' F°ll° w i n g the same analysis as for asters, 
the instability of vortices occurs for 

v ^ m \ 2 4r? + 7 1 (z/ 1 + l) 2 



The resulting stability diagram in the (C^A*, 5i^)-plane is displayed on Fig. 2. We have 
assumed on this diagram that V\>\. There are three regions on this diagram, a region 
with positive 5K and positive (Afi where the asters are stable, a region with negative 5K 
and positive (Afi where vortices are stable and a region with large negative (Afi where 
spiral defects are stable. Along the line 5K = 0, spirals are always stable as studied in the 
previous paragraph. We do not give here a general description of the spiral defects, as in the 
case where SK = 0, we expect that there is a dynamic selection of the orientational angle 
and that the defect is rotating. 
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D. Effects of friction with the substrate 

So far we have considered the point defect in a real two-dimensional space. A real active 
gel interacts with both the surrounding solvent and the substrate. We do not discuss the 
hydrodynamic interactions between the gel and the solvent here. We now consider a thin 
quasi-two dimensional gel layer and introduce a friction force at the substrate proportional 
to the local velocity. Taking into account the rotational symmetry, the friction force modifies 
the ortho-radial (9) component of the force balance that becomes 

— h -a rg = £v e , (69) 

or r 

This equation, which is a balance between the friction force and the tangential stress can 
be integrated for a rotating spiral defect, yielding the following expression for the tangential 
stress as a function of the velocity. 

(t,. h = X C dr' r' 2 v 9 (r') (70) 



r 2 Jo 
In the absence of elastic constant anisotropy, SK = 0, we can derive the velocity profile 

vg(r) in an infinite system from Eq. (|70|) as follows: We look for a solution of Eq. (|52|) with 

h± — and ip =const.. As in the absence of friction, the polarization angle can have one 

of the 4 possible values satisfying cos2^ = \jv\. The velocity field is then obtained from 

Eq. (|55jl that can be transformed into the following ordinary differential equation for vq: 

d 3 / v e ldv e \ 2 . \ i^\ 

where z = r/Xf and the the friction length is given by A/ = (4r/ + jxUi sin 2 2ipo) 1 ' 2 /(2£ 1 ' 2 ). 
The velocity is then given by 

v $ {r) = 2 tao X f (K 1 {y)-^Y (72) 

where K\ is the modified Bessel function of the second kind defined in Ref. [3^ 

At short distances (r <ti Xf), the dissipation is dominated by the viscosity rj and the 
velocity field is the same as in the absence of friction. At large distances, the dissipation is 
dominated by the friction on the substrate £ and the velocity vg decays to zero. The stability 
diagram of the defect in the presence of substrate friction is very similar to the diagram of 
Fig. (3) where the finite size R would be replaced by the friction length Xf. 
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VI. DISCUSSION 

We have introduced in this manuscript equations which describe the long wavelength and 
low frequency behavior of active gels. Although we have written the equations specifically in 
the case where the activity is due to motor proteins, they should apply, in their principles, 
to all gels in which a permanent source of dissipation is at work. Such gels define a new class 
of materials. For instance, a conventional "physical" polar gel, absorbing a high frequency 
ultrasonic wave should obey, in the low frequency, long wavelength limit, the set of equations 
proposed here. Our main motivation however, is the construction of a generic theory for 
characterizing quantitatively the properties of Eukaryotic biological gels. One could object 
that generalized hydrodynamic theories involve a large number of parameters and are thus 
not very useful. Their merit is to involve the smallest number of parameters required for 
a relevant description of the systems under consideration and to describe in a unified way 
all long length scale and long time scale situations which are otherwise seemingly uncon- 
nected. The nature of most parameters is already well-known for gel or polymer rheology 
or for liquid crystal physics, and their measurement techniques can straightforwardly be 
transfered to active gels. This is transparent for translational, orientational elastic moduli 
and viscosities. Coefficients linking shear flow and polar orientation are less familiar to the 
general public but are well known to liquid crystal physicists, and their measurements are 
not difficult a priori. There are six bulk and two surface additional parameters, compared 
to a passive polar gel. Among the bulk quantities one is equivalent to the motor velocity 
on the actin filaments (A), and another one is the depolymerization rate kj of the gel in the 
bulk. Both quantities are directly accessible to experiment. Surface polymerization terms 
can be extracted independently from biochemical data and although they generate a new 
interesting physics (|3ll34|). they do not introduce uncertainties in the description. 

Among the four remaining parameters, three ((, ( and (' bear essentially the same physics 
(i.e: the activity implies either spontaneous motion or spontaneous stress), and the last one 
Ai measures the active rate of change of polarization. Hence, there are only two relevant 
additional parameter, namely Ai and ( with respect to a passive polar gel. We show, in the 
two examples developed here, that they change profoundly the behavior of gels: behaviors 
which in the absence of energy input would be static, become dynamic. The structure of 
the relaxation modes which would be entirely over-damped in the long wavelength limit, 
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now can support propagative waves 



24l.l3l|. spiral disclinations rotate permanently. Many 
other consequences have to be unraveled. Knowing that it has taken more than ten years to 
investigate the properties of generalized hydrodynamical equations relevant to liquid crystals, 
it is likely that a similar number of years will be necessary in this case as well. Our bet is 
that it will help us understand the complex behavior of the slow dynamics of the eukaryotic 
cytoskeleton in a robust way, not depending on the details of the involved proteins. 

Our theory can be extended in several ways. First we have ignored the permeation process 
of the bulk fluid through the gel. This is legitimate in the long time limit in most geometries 
but not all. The inclusion of permeation in our equations would be straightforward. We have 
also used the simplest visco-elastic gel theory: experiments performed on cells suggest the 

existence of scale invariant visco-elastic behavior which could also be included in the theory 

"I 

35| . At last we have assumed that the driving force was " small" and constant both in space 

and time. The spatial and temporal invariance make sense in a cell since the ATP produc- 
tion centers seem to be abundant and evenly distributed. In vitro experiments may mimic 
these conditions on a limited time scale, yet large compared to most phenomena of interest. 
Furthermore there is no a priori difficulty in introducing ATP production and consumption 
in the equations. More severe is the last limitation: our equations are valid for A/i small 
compared to thermal energy, whereas it is of the order of ten time that value in real life. 
The extension of our theory to large A/i is a real challenge and will probably require some 
feedback from experiments. A brute force expansion in higher powers of forces and fluxes is 
just totally impractical and would leave out the known and probably leading exponential de- 
pendences of polymerization/depolymerization rates on stress. Other challenging questions 
concern the nature of noise which will pick up strong out of equilibrium contributions, and 
possible dynamical transitions for instance in the motor behavior 36]. For these reasons, we 
think that is it wise to start with the well controlled approach that we propose here and to 
develop the corresponding experiments. This will help us familiarize with this new physics 
before getting to regions of phase space were we lack guiding principles. 

An important application of active gels is the problem of cell locomotion on a substrate 
J37|. In a forthcoming publication, we will discuss how the interplay of polymerization and 
activity can induce the motion of a thin layer of an active gel on a solid surface. The 
physics of this gel motion coupled to the properties of bio-adhesion provides a basis for an 
understanding of the locomotion of cells such as keratocytes [38] . 
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APPENDIX A: CYLINDRICAL COORDINATES 

In polar coordinates, the rate of strain tensor for a rotationally symmetric velocity profile 
is given by 

d (Al) 

(A2) 

(A3) 
z ar v t / 

Furthermore, 

1 d 

(A4) 

If the system is invariant by rotation, the local force balance in cylindrical coordinates is 
written as 

;!(«#)-;* = £ («) 

hi ( rV ") = ° • (A6) 

In an incompressible system, the pressure II is a Lagrange multiplier necessary to impose 
the incompressibility constraint. 

APPENDIX B: VISCOELASTIC POLAR GEL 

The form of Equations ()17jl. (jTSjl . (|23jl and PB} can be understood by discussing the 
general dynamics of the polarization field p a in frequency space. As a function of frequency, 
we write to linear order 

h a I 

iup a = — + XxPaAfl - — — {v X ppU a p + Vip a Upp)... (Bl) 

71 1 + IUOT 

where the tilde denotes a Fourier amplitude at frequency uj. The first two terms of the right 
hand side correspond to a straightforward linear response theory of a ferroelectric liquid 
crystal. Similarly, the last term provides for u = the reactive linear coupling term allowed 
by symmetry. For finite frequency furthermore, this term takes into account that the gel is 
elastic on short time scales or large uj. This can be seen as follows: in the case of an elastic 
gel, the free energy F of the polarization field contains elastic terms of the form 

d d xp a e a ppp (B2) 



2r 
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where e a p is the strain tensor characterizing deformations of the gel. This leads to a mod- 
ification of h a by h e ^ = —6F e i/5p a - Since in frequency space, u a p = ioje a p, Eq. (jBljl does 
describe the correct elastic limit for large u Separating Eq. (|B1J) in reactive and dissipative 
parts, using the opposite signatures with respect to time reversal (or equivalently uj —* —uj) 
allows us to identify P^ and P T a . Using 1/(1 + iur) = 1 — iur/(l + uj 2 t 2 ) and replacing iu 
by D/Dt, we can write the reactive and dissipative parts of P a as given by Eqns. (fTHjl and 
(J25J) . Symmetry relations now impose the correct equations for a^ and a r a/3 as given by 
Eqns. (JT7J) and 
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a) \ 1 / 





FIG. 1: Schematic representation of the geometry of point defects of topological charge one in 
two dimensions. Displayed are the orientations of polarization vectors as arrows for (a) aster, (b) 
vortex and (c) spiral. 
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FIG. 2: Stability diagram of asters, vortices and spirals which are topological defects in an active 
gel of polar filaments. Asters are stable for 5K > in the region with actively generated stresses 
£A/z larger than a critical value C^-Hc ■ This critical value is negative, corresponding to contractile 
stresses in the gel. Vortices are stable for 5K < and £A/i > C^-fi^ . For other parameter 
values, rotating spirals occur via a symmetry breaking dynamic instability. Here, 5K = 8K/K is a 
dimensionless ratio of two elastic moduli and £A/z = -R 2 (A/x/-K" a dimensionless measure of active 
stresses. Note that both rotation senses occur with equal probability on symmetry grounds. The 
diagram was evaluated for the choice 77/71 = 1 and v\ = 2 of Onsager coefficients of the system. 
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